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Abstract 

"Relevation"  is  the  name  given  by  Krakowski -flST^)  (Rev.  Franqaiee 
Autonzt.  Inform.  Recherche  Op.,  7,  Ser.  V-2,  108-lZtt^ to  the  distribution  of 
failure  time  of  a  replacement  from  an  aging  stock.  The  present  authors 
(Johnson  and  Kotz  (49?9t''T98T)  (WEE  Trane.  Reliability,  R-28 ,  292-299; 

American  Journal  of  Mathematical  and  Management  Sciencei ,  If have  extended 
this  concept  to  include  (i)  hierarchal  replacement  systems  and  (ii)  dependence 
between  lifetimes  of  original  and  replacement  items.  In  this  paper,  we  present 
some  further  developments,  including  first  steps  towards  a  synthesis  of  (i)  and 

(ii).  . 

Key  Words  and  Birases:  Relevation;  Dependence;  Aging  in  Stock;  Aging  in 

Service;  Farlie-Gunbel-Morgenstem  Distribution; 
Hierarchal  Systems. 
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1.  Introduction. 

The  concept  of  a  "relevation"  was  introduced  by  Krakowski  (1973)  to 
correspond  to  replacement  from  aging  stocks.  Suppose  that  when  an  element  with 
survival  distribution  function  (SDF)  S^(t)  fails,  it  is  replaced  by  an  element 
with  SDF  S2 Ct)  taken  from  stock.  If  failure  occurs  at  time  tj,  then  the 
conditional  SDF  of  the  replacement  element  is 

S2 | 1 ^t I tj)  *  S2(t)/S2(t1)  , 

and  the  SDF  of  time  of  failure  (T)  of  the  replacement  element  is 

St(T)  «  S^T)  ♦  S2(t)  Jj  f1(t){S2(t)}'1  dt  ,  (1) 

where  fj(t)  *  -dSj(t)/dt  is  the  probability  density  function  (PDF) 
corresponding  to  Sj (t) . 

Sp(r)  is  called  the  relevation  of  Sj(t)  with  S2(t). 

This  concept  has  been  extended  to  situations  wherein  there  is  a  hierarchal 
system  of  replacement,  and  wherein  the  lifetimes  of  elements  are  not 
independent  (Johnson  and  Kotz  (1979,  1981)).  In  the  present  paper,  we  study 
combinations  of  extensions  of  this  kind,  and  also  give  exanples  in  which 
differential  rates  of  aging  in  storage  and  in  service  are  allowed  for. 

2.  Dependent  Relevations  in  a  Hierarchal  System. 

Consider  a  set  of  n  elements  with  a  simple  replacement  element  £2. 

This  might  occur,  for  example,  when  £2  is  called  into  service  as  a  back -13) 
only  when  . there  is  a  failure  among  {£j}.  This  system  (conprised  of  (£j)  and 
C2)  would  quite  possibly  be  only  a  small  part  of  a  conplete  hierarchal  system; 
study  of  such  relatively  simple  systems  is  an  essential  prerequisite  to  study 
of  more  elaborate  systems  of  which  they  could  be  a  part. 


We  denote  the  failure  times  of  the  elements  of  {5^}  by  Tjj,Tj2»***»Tjn» 
and  of  ^2  by  T2"  j°*nt  surviva*  distribution  (SDF)  of  these  (n  +  1) 


failure  times  is 


SJ2^~l;t2^  *  P(j2iCTlj>tlj)  0  ^2>t2^ 


(tj  *  (tjj . tln)).  Joint  SDF  oF  Ti  =  n'in(T11,...,Tln) --the  failure 

time  of  the  first  member  of  {E^}  to  fail  --  is 

- - 1 

Sl2(ti'»t2^  =  S12^tl~*  *  t2^  rNTrs  WA&i 

r"-'U:  TAB  r  / 


where  1’  =  (1,1,...,  1);  and  the  SDF  of  is 


S^tj)  *  S12(tj  1’ ;  0)  . 


The  SDF  of  the  replacement  element  is  then 


•  -;nced 
<hjJtifieatio.»> 


rT  ^2 1 1  (T  1 1  ^  ^r^J- — 


where  fj(t)  *  -dS1(t)/dt  and  S^ 1 1 (T 1 1)  “  PC^  >t |T^*t) . 

In  general,  explicit  evaluation  of  (5)  leads  to  rather  complex  formulas. 
As  a  compromise  between  realism  and  simplicity,  we  will  (similarly  as  in 
Johnson  and  Kotz  (1981))  use  a  generalized  Farlie-Gumbel-Morgenstem 
distribution  (Johnson  and  Kotz  (1975))  as  a  special  case  for  detailed 
examination.  We  suppose 
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n  n 


♦  8  [  ?  (1  *  S(t  .)}{1  -  SCt,.,)}]  .  (6) 

j  <  j  *  13 


The  parameters  a  and  e  have  to  satisfy  certain  conditions  for  the  distribution 
to  be  a  proper  SDF.  These  conditions  can  be  sumnarized  as 


|(2r-n)a|  -  >s  { f  2r -n) 2  -  n)6  si  (r  *  0,1,.. . ,  n)  , 


or  equivalently. 


(7) 


r'Jo)  -  h(r'Z-n)B  <  1 


r'  *  0,2,...,n  for  n  even 
r*  =  l,3,...,n  for  n  odd 


The  marginal  SDF  of  each  T^  is  S(t) ;  the  marginal  SDF  of  T^  is  S^Ct).  The 
joint  SDF  of  ^2 !»***»  "^in  *s 


St  (t,)  =  {  n  S(t, .) }  [1  ♦  6  Y  l  {1  -  S(t,.)Hl  -  set,.,)}]  ,  C8) 

il  1  j=1  j  <  j  *  lJ  13 

which  is  an  exchangeable  distribution.  This  is  likely  to  be  a  reasonable 
assunption  if  there  are  no  systematic  differences  among  the  conditions  to  which 
the  n  elements  of  (Sj)  are  exposed. 

The  SDF  of  Tj  (time  to  first  failure  in  {£j})  is 

Sj(tj)  *  (SCtpf  tl  ♦  JsnCn-l)B{l  -  Sftj)}2]  ,  (9) 

and  the  corresponding  density  function  (PDF)  is 
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,n-l 


fl<V  “  "  [1  +  %(n  -l)8(n  -  2(n  +  l)S(tj)  +  (n+2)  [SCtj)]2}] 


-  nfCtjHSCtj)}""1  [1  ♦  Mn-DBd-SCtj))^  -  (n  +  2)S(t1)}]  .  (10) 


The  joint  SDF  of  Tj  and  T2  is 


S12(tl’t2)  *  [1  *  na{l  -  SftjjH  1  -  S2(t2)} 

+  %n(n-l)B{l  -  S(tp}2]  , 


(11.1) 


with  joint  PDF 


f12(tl’V  *  nf1(t1)f2(t2){S(t1)}n'1  [1+  a(n  -  (n+DSCt^Hl  -  2S2(t2)} 


+  Mn-1)B(1  -  SCtjjHn  -  (n+2)S(t1)}]  . 


(11.2) 


The  conditional  PDF  of  T2>  given  Tj,  is 


f  ft  It  1  -  f12('tl,t2-) 

2|1(  2(tl}  ~  -  fjCtj) 


f2(t2) 


l  +  a(n-  (n+IJSftpHi- 2S2(t2)}  +  >s(n-l)B{l  -  Sftj)}  (n  -  (n^JSftj)} 
1*  *s(n-l)B{l  -  Sftj)}  (n  -  (n^jSftj)} 


(12.1) 


and  the  conditional  SDF  is 


c  r*  I*  ]  -  5  ft  l  +  «(n-  (n+l)S(t1)Hl-S2(t2)}+Jj(n-l)B{l-S(t1)}{n-(n+2)SCt1)} 
2|1  2'  1  "  2Ct2  l*%(n-l)BU  -  SCt^Hn  -  (n+2)SCt1)} - 


(12.2) 
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Hence  from  (5)  the  SDF  of  time  to  failure  (T)  of  the  £2  (replacement)  element 
is 

Sj-Ct)  -  {S(r)}n  (1  ♦  Jjn(n -1)6(1  -  S(x)}2] 

♦  n  /J  f(t)(S(t)}n'1  [1  ♦  *s(n  -1)6(1  -  S(t)Hn  -  (n+2)S(t)}] 

S2(t) 

X  — - 

s2(t) 

1  ♦  all  -  S2(t) }  (n  -  (n+DSCtj) }  ♦  %(n-l)6(l  -  Sft^} (n  -  (n+2)S(tj)} 
1  ♦  a(l  -  S2(t)}  (n  -  (n+DSCtj)}  +  %(n-l)B(l  -  S(t)Hn  -  (n+2)S(t)} 

(13.1) 

As  a  further  specialization,  we  will  take  S2(t)  «  S(t).  Then  making  the 
substitution  s  =  S(t)  in  the  integral  leads  to 

Sj.(t)  *  (S(T))n  [1  ♦  %n(n -1)6(1  -  S(t))2] 

♦  nS(r)  /*(T)  sn'2[l  +  Js(n- 1)6(1  -s)(n  -  (n+2)s» 

l  +  ct(l  -  S(t))  (n  -  (n+l)s)*>i(n-l)6(l-s)(n  -  (n+2)s)  . 

1  ♦  a(l-s)(n  -  (n+l)s)  +  Jj(n-l)6(l-s){n  -  (n+2)s) 

(13.2) 

Note  that  £j.(t)  is  always  the  same  function  of  S(t),  whatever  the  functional 
form  of  the  latter. 

The  integral  in  (13.2)  can  be  evaluated  in  terms  of  elementary  functions 
(see  Appendix) ,  though  the  expression  is,  in  general,  rather  complicated.  For 
the  special  case  o  ■  0  (so  that  replacement  lifetime  is  independent  of  the 
lifetimes  of  the  initial  elements) , 


dt 
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Sj-Or)  -  ($(T)}n  [1  ♦  »sn(n-l)e{l  -  S(t)}2] 

♦  nS(x)  /*(t)  sn"2[l  ♦  %(n  - 1)8(1  -s){n  -  (n+2)s)]<Js 

«  S(t)  [n  -  {SCO)0'1] 
n-1 

+  Mli[l  -  MSCO)"*1  [n(n  +  l)  -  2(n2 -l)S(t)  +' n(n-l){S(r)}2]] 
n+1 

■  \  *  V  •  <14> 

Some  numerical  values  of  A  and  B  are  shown  in  Table  1. 

n  n 

TABLE  1 

Values  of  and  B^  for  S^.(t)  *  A^  +  Bn  8  when  a  *  0. 


n  2  3  4  5 


S(T)\  « 

Bn 

A 

n 

B 

n 

Bn 

An 

B 

n 

0.2 

0.360 

0.0341 

0.296 

0.0410  | 

0.266 

0.0377 

0.250 

0.0328 

0.4 

0.640 

0.0288 

0.568 

0.0475 

0.524 

0.0S46 

0.497 

0.0547 

0.6 

0.840 

0.0128 

0.792 

0.0269 

0.757 

0.0381 

0.731 

0.0487 

0.8 

0.960 

0.0021 

0.944 

0.0054 

0.930 

0.0094 

0.918 

0.0132 

(NOTE:  Possible  values  of  8  are  limited  by  -2n  *(n-l)  *  s  6  s  2n  *.) 

If  a  «  0,  the  value  of  St(t)  exceeds  that  given  in  (14)  by 

naSfTl  f1  sn~2n*  %(n-l) 8(1-5) (n-  (n+2)s)Hs  -  S(t)  )  (n  -  Cn+l)s)  . 

J  jS(t)  l  ♦  a(l-s)(n  -  (n+l)s)  +  h(n-l)B(l-s)  {n  -  (n+2)s)  05 

If  S(t)  >  n/(n*l),  this  will  be  of  opposite  sign  to  a  since 

{s  -  S(x)Hn  -  (n+l)s)  <  0  for  S(x)  <  s  <  1  . 


This  gives  some  qualitative  appreciation  of  the  effect  of  non-zero  a. 

To  see  how  these  results  can  be  used  in  dealing  with  more  conplicated 
systems,  consider  the  subsystem  set  out  schematically  in  Figure  1.  The 


FIGURE  1 

component  £3  is  related  to  the  two  £2's  ^  tlie  same  W  35  each  ^2  is 
related  to  the  two  Sj’s  belonging  to  it.  £3  comes  into  operation  at  the 
first  failure  of  the  two  £2's.  The  relevation  SDF  Sy(x)  given  by  (13.1)  now 
plays  the  role  previously  played  by  the  SDF  of  each  (S(t)).  Of  course,  if 
the  resulting  formula  for  the  SDF  of  the  £3  replacement  were  to  be  written  out 
explicitly,  it  would  be  very  lengthy  indeed.  However,  the  computation  of 
nunerical  values  is  quite  straightforward.  It  can  be  carried  out  in  stages: 
first  finding  £j.(x)  and  then  combining  it  with  the  SDF--Sj(t),  say --of  the 
£3  components.  (The  case  of  independent  conponents  has  been  treated  by 
Johnson  and  Kotz  (1979).) 


3.  Differential  Rates  of  Aging  in  Storage  and  in  Service. 


It  is  likely  that  an  element  will  be  under  less  stress  while  in  storage 
than  in  service.  This  can  be  allowed  for  by  regarding  a  period  t  in  storage  as 
equivalent  to  a  period  g(t)  in  service  (usually  with  g(t)  <  t).  Equation  (1) 
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is  then  replaced  by 


Sj(t) 


S^t) 


s2(g(t)+  T  '  t) 
S2(g(t)) 


dt  , 


(15) 


and  (2)  by 


Sj(T) 


fo  flW 


S2|l(g(t)  +  T-t|t) 

s2)1(g(t)|t) 


dt 


(16) 


One  might  reasonably  take  g(t)  *  yt  —  so  that  1  year  in  storage  is 
equivalent  to  y  years  in  service.  Usually  one  would  take  y  <  1,  but  this  is 
not  essential.  Even  if  it  is  assumed  that  Sj(t)  *  S2(t)  *  S(t),  it  is  no 
longer  true  that  depends  only  on  S(t)  .  As  an  example,  we  take  y  <  1 

and 


S(t)  =  1  -  Kt  (Ostsrs  k ’*)  . 
Then  from  (15) , 

¥*>  ■ 1  •  «  *  C « * 

*  1  -  tcry  *  -  (1  -  Kry)y 


log(l  -  tcry)  . 


(17) 
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Some  numerical  values  are  given  in  Table  2. 


TABLE  2 

Values  of  St(t)  when  S(x)  »  1  -  ict  has  specified  values 
(y  *  service  years  equivalent  to  one  year  in  storage) . 


smV 

0.2 

0.4 

0.6 

0.8 

1.0 

0.2 

0.661 

0.639 

0.611 

0.562(5) 

0.522 

0.4 

0.812 

0.804 

0.794 

0.781 

0.767 

0.6 

0.918 

0.915 

0.913 

0.910 

0.906 

0.8 

0.980 

0.979 

0.979 

0.979 

0.979 

(The  last  column  corresponds  to  no  difference  between  aging  in  storage  and  in 
service.)  The  effect  of  differential  aging  is  not  very  great,  especially  at 
the  longer  durations. 
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The  integrand  in  (13.2)  can  be  expressed  in  the  form 


[a  {(s-b  )2  ♦  c  l]*1  T  d  .sn+J'2 
n  n  n  n,] 


Using  the  change  of  variable  y  «  s  *  b  ,  the  integral  becomes 


“n1  Ism- b  {,L  dn.i(y*b7.,"‘2*:'K>,Z*(:n)'1  d>' 


n  j=0 


n,j(y*br.)"'2,:iKy2‘Cn>'1  <* 


_  -1  r  j  n‘?^  rn-2+j..  n-2+j-h  f1-bn  Ji,  2  .-1  , 

'a"  Jo  "•i  Jo  ‘  »  )b"  !SM-bnyb  *c„>  ^ 

-*;lb"'2  j0  dn,jbn  X  r^)b;b[Ih(l-bn)  -  IhCS(T)-bn)i  , 


where 


Xh(y)  =  ft  yh(y2’,‘cn)  '1  dy 

(h  odd) 

(h  even,  cn  >  0) 

(h  even,  c^  <  0)  . 

[For  exanple,  Iq(y)  b  l/*1^  tan’1  y/Jc^  if  cfi  >  0; 

[1/(2/^ )] log C(y-/^)/(y+/^))  if  cn  <  0;  IjCy)  =  h  log(y2+cn); 

I2(y)  “  y  -  cnI0(y);  Ij(y)  -  >sy2  -  c^Cy);  I4(y)  -  (l/3)y3  -  cnI2(y);  and  so 
on.) 


l°g(y2  +  cn) 

[(h-3)/2)  (-c  )*  .  h  .  -  n  , 

j0 

*>. 

7  n 
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The  constants  a  ,  b  ,  and  c  are : 

n  n  n 

an  =  (n  +  l)a  +  *s(n  - 1)  (n  +  2)B  ; 
bn  =  *s{(2n  +  l)a  +  (n^-ljBia’1  ; 
cn  =  tan  '  l/4(a  +  (n-l)6}2]a^2  . 

(Note  that  if  a  =  0,  an*  ^(n -1)  (n  +  2)6;  =  (n  +  l)/(n  +  2); 

cn  «=  2{(n-l)(n  +  2)B}'1  -  (n  +  2)'2.) 

The  egressions  for  the  coefficients  d^  ^  are: 

dn  0  =  {1  +  !sn(n-l)B)[l  +  nail  -  S(x)}  +  *sn(n-l)B]  ; 

d  ,  *  -2(n2-l){l  +  ^ n(n-l)B)B  -  (n  +  l)a{l  +  (3/2)n(n-l)6}{l  -  Sfx) } 
n ,  i 

dn  ,  =  (n  -l)(n  +  2)B  +  %(n -l)(3n2  +  6n  +  2)  [a(l  -  S(r)}  +  (n-l)B]B  ; 

d  _  =  -(n2  -1)  (n  +  2)  [*sa{l  -  S(x)}  +  (n-l)B]B  ; 
n,j 

d  .  =  (l/4)(n-l)2  (n  +  2)2  B2  . 
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"Relevation"  is  the  name  given  by  Krakowski  (1973)  {Rev.  Fran^oise 
Automat.  Infom.  Recherche  Op. ,  7 ,  Ser.  V-2,  108-120)  to  the  distribution  of 
failure  time  of  a  replacement  from  an  aging  stock.  The  present  authors 
(Johnson  and  Kotz  (1979;  1981)  (IEEE  Trane.  Reliability,  R-28,  292-299; 
American  Journal  of  Mathematical  and  Management  Sciences,  7(2)))  have  extended 
this  concept  to  include  (i)  hierarchal  replacement  systems  and  (ii)  dependence 
between  lifetimes  of  original  and  replacement  item*..  In  this  paper,  we  presen 
some  further  developments,  including  first  steps  towards  a  synthesis  of  (i)  an 
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